Phase-space correlations of chaotic eigenstates 
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It is shown that the Husimi representations of chaotic eigenstates are strongly correlated along 
classical trajectories. These correlations extend across the whole system size and, unlike the corre- 
sponding eigenfunction correlations in configuration space, they persist in the semiclassical limit. A 
quantitative theory is developed on the basis of Gaussian wavepacket dynamics and random-matrix 
arguments. The role of symmetries is discussed for the example of time-reversal invariance. 
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Chaotic eigenfunctions and in particular their localiza- 
tion and correlation statistics are a topic of continuing in- 
terest [J 0> 0, 0> El • Applications include classical, 
mesoscopic and pure quantum systems such as optical, 
mechanical and microwave resonators 0, , electron 
transmission and interaction in chaotic quantum dots 
[HI , and decay and fluctuations of heavy nuclei . 
One foundation of eigenfunction statistics is the random- 
wave model of Berry 1] which is essentially equivalent 
to random-matrix theory (RMT) 2]. Within RMT the 
eigenfunction components in an arbitrary basis are un- 
corrected Gaussians. Current research is frequently aim- 
ing at deviations from RMT due to the specific dynamics. 
Prominent examples are scarring by periodic orbits p| or 
long-range correlations 6]. As classical dynamics takes 
place in phase space, representations of eigenstates via 
Husimi or Wigner functions seem appropriate, and re- 
cently some of their statistical properties have attracted 
a lot of attention 7]. Surprisingly this does not apply 
to dynamically induced correlations although numerous 
studies of the corresponding spatial correlations demon- 
strate their relevance 0, U E]> an d although there are 
systems where a direct relation between phase-space cor- 
relations and measurable quantities must be expected. 
For example, in optical resonators 8] the power emit- 
ted at a certain point of the boundary depends strongly 
on the angle of incidence of the wave (total internal re- 
flection). Also in quantum dots effects of eigenfunction 
directionality have been measured using tilted leads [l2| . 

In this paper we analyze for the first time dynamical 
correlations between points in phase space. Our results 
are surprising in view of the fact that spatial two-point 
correlators of eigenfunction amplitudes or densities van- 
ish in the semiclassical limit Ti — > for any x =^ x' . In 
contrast we find strong and semiclassically persistent cor- 
relations between phase-space points £ ^ £' (£ = (x,p)) 
where the distance between x and x' can be of the or- 
der of the system size. This is no contradiction since 
(i/j(x)il)(x')) = does not imply that ip(x) and ip(x') are 
statistically independent; it just means that there are no 
linear correlations in configuration space. In other words, 
although the existence of dynamically induced correla- 
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FIG. 1: (a) An eigenstate of a chaotic torus map quantized 
with N = 256 is shown in Husimi representation. A classical 
trajectory £t = (xt,pt) of the same map is marked by circles 
with label t. Inset of (b): In a scatter plot the values of the 
Husimi density at £o and £i cluster near the diagonal which 
indicates strong correlation, (b) The correlation coefficient is 
shown with circles as function of N. The dashed line is the 
prediction of Eq. (|13f> . The squares show that spatial density 
correlations between the corresponding positions xo and x\ 
vanish as N — > oo. 



tions cannot be a matter of the chosen basis they turn 
out to be most relevant in phase space. 

In our approach we make use of methods which proved 
successful in studies of eigenfunction scarring. We will 
see that Gaussian wave packet dynamics |l5| can also be 
applied to study eigenfunction properties which are not 
at all related to periodic orbits. As in nonlinear scarring 
|sj we have to supplement short-time dynamics by RMT 
in order to get a complete theory. However, the new 
context requires to do so in a technically different way. 

We introduce our main result with a numerical exam- 
ple. Fig.^Ji shows the Husimi representation of an eigen- 
state of a quantum kicked rotor, 



U 



e -27riV(x)/h e -2iriT(p)/h 



(1) 



The unitary matrix U has dimension N = h 1 and quan- 
tizes a classical map on the torus (x + 1 = x, p + 1 = p), 



x t +i x t 



T\p t ) 



Pt+i Pt 



V\x t+1 ). (2) 
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We choose T(p) = p 2 /2 and a potential for which this 
map is globally chaotic [l6|. Therefore the eigenstates 
of U fill the entire phase space more or less uniformly. 
The fluctuations on this background are the object of 
our interest. In Fig. we observe patches of high or 
low density which are typically of the size of a Planck 
cell h. Besides this well-known local correlation 7], the 
figure suggests the existence of long-range correlations 
along the trajectories of the underlying classical dynam- 
ics. One such trajectory £ t (not periodic!) is indicated 
by circles. The point £o is close to a maximum of the 
density, and we observe that also at the points the 
Husimi density has values above the average, at least for 
— 1 < t < 2. A scatter plot of the densities at £o and £i 
shows that this is no coincidence (inset of Fig.GJ)). The 
corresponding correlation coefficient (circles in Fig. 
is as high as 0.94 although the points £o and £i are far 
apart and not related by any symmetry. This is in strong 
contrast to the absence of long-range correlations in po- 
sition or momentum representation in the semiclassical 
limit, which is expected from previous studies 0, H, 
and confirmed in Fig.EJ} f° r °ur model (squares). 

The existence of dynamical correlations in phase s pac e 
can be derived from Gaussian wave packet dynamics 15]. 
For the special case of the quantum map Eq. (QJ this 
standard approach is formulated as follows. The wave 
function of a general Gaussian state centered in phase 
space at £o = (xo,Po) has the form 

0(fo, A£, AT, s- x) = j\f e - A ^A x - x o) 2 +ivo{x-x Q )/n ^ 

where A, M are complex. The prefactor M accounts for 
normalization and an overall phase while 



Aaz,s = (2A*r l + s(2i/ft) 



A^ 



h 2 /4 (4) 



depends on the uncertainties of position and momentum, 
A£ = [A x ,A p ], and on a sign s = ±1. We assume that 
both, A x and A w , are much smaller than any relevant 
classical scale 18]. In this case and within the semiclas- 
sical approximation (stationary phase approximation for 
N — > oo ) application of the propagator preserves the 
Gaussian form of the wave function, 

U\g(£o,&Zo,Mo,s )) = |s(£i,A£i,M,si)>. (5) 

In this equation £i is the classical iterate of £o- Position 
and momentum uncertainties transform as 



A 2 xl = A 2 - 2a r 'JAl A* - h 2 /4 + A 2 pQ Tf (6) 
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and also for A/*i and si semiclassical expressions are 
found. In a chaotic system the growth of the uncertain- 
ties under repeated application of Eqs. (0, J7J) will be 
dictated by the classical Lyapunov exponent A after a 



short initial period because then h 2 /4 <C A xt A^ can be 
neglected and one is left with equations that are essen- 
tially classical. Hence, for A^ ,A^ Q ~ Ti the uncertain- 
ties grow to values ~ 1 after a time of the order of the 
Ehrenfest time £e ~ A -1 Info and then the approximate 
Gaussian wave packet dynamics breaks down. 

We can now project the semiclassical iterates \gt) = 
\g(£t, A^t, A/"t, s t )} of a Gaussian \go) onto an eigenstate 
\n) of the propagator, U \n) = exp(ie n /h) |n), and obtain 



(n\g t ) = exp(i£ n t/h) (n\g ) 



(8) 



for times up to the Ehrenfest time. Already from this 
simple equation it is obvious that the projections of eigen- 
states into phase space along a classical trajectory can- 
not be independent of each other and that there must 
exist cross-correlations between eigenvalues and eigen- 
states. Note that Eq. JHJ) is a semiclassical identity which 
has to be satisfied by every single eigenstate even with- 
out averaging. However, the relevance of this equation is 
diminished by the fact that in general \go) and \g t ) are 
Gaussian wave packets with completely different widths 
and shapes. Therefore it is not immediately possible 
to extract meaningful correlation functions relating, e.g., 
Husimi densities at different points in phase space. This 
problem will be addressed in the following. 

For the sake of simplicity we will consider coherent 
states with equal position and momentum uncertainties, 



%W = (ft* - )' 



-1/4 -{x-xqY /2h+ip (x-x )/h 



(9) 



The Husimi amplitude of a state \ip) at a point £ is 
^(0 = ( a fl^) an d the Husimi density is H^(£) = 
1^(0 1 2 - We are mainly interested in the correlation co- 
efficient (normalized covariance) for Husimi densities at 
two different points £ and It is given by 



(SH n (QSH n (?)) n 
^/(5 2 H n (0) n (S 2 H n ^)) r 



(10) 



where (-) n = ^V _1 X^=i(') denotes the average over 
all eigenstates \n) and 5H n (£) = H n (£) — N~ x is the 
deviation of the Husimi density from its mean value 
(iJ n (£)) n = TV -1 . It will be instructive to consider in 
parallel also a suitably defined time-dependent correla- 
tor of Husimi amplitudes 



c t 



(e , ;f) = l^(e +te " t/ ^;(O^(0>r 



(11) 



which involves the eigenphase e n . Fig. illustrates the 
main properties of the correlation functions ([TUj). ([TT]) 
which we shall explain below semiclassically. Fig. 
shows for t = 1 that ct(£o;0 as a function of £ is con- 
centrated at the classical iterate The magnitude of 
the correlation between these two points depends on the 
classical trajectory. This dependence is shown in Fig. [2^ 
by plotting ci(£o;£i) as a function of £o- Figs.Et>,d con- 
tain the same information for the correlator Ch whose 



3 




FIG. 2: (a) For fixed £o the correlator of Husimi amplitudes 
ci(£o;£) i s displayed (N = 1024). The point £o and its classi- 
cal iterate £i are marked by circles, (b) Same for the density 
correlator Ch(£o', 0- Labeled circles show the classical trajec- 
tory £t- Unlabeled circles show the time-reversed trajectory 
£t obtained by the transformation flS) . (c) ci(£o; £i) is shown 
as a function of £o- The magnitude of correlation is approxi- 
mately constant along lines x\ — xo + po = const, (d) Same 
for Ch(£o;£i)- On top of a background depending on x\ one 
observes smooth fluctuations and relatively sharp lines of ex- 
cess correlation which correspond to Eq. (J2UJ) . 



structure is more complicated. In Fig. we see that the 
Husimi density along a classical trajectory £ t is correlated 
over a short time, as expected from Fig. ^ I n addition 
we observe also a strong correlation to the time-reversed 
trajectory and some fluctuating background of weakly 
positive or negative correlations covering the entire phase 
space. The dependence of Ch (£o?£i) on the position in 
phase space (Fig. 01) resembles the characteristic pat- 
tern observed already for the correlator of amplitudes in 
Fig. I2t- However, it is blurred by fluctuations and has 
in addition sharp lines of particularly high correlation 
which we shall relate to the time-reversal symmetry of 
our model. 

It is straightforward to obtain a semiclassical approxi- 
mation to Eq. ((11)). Using the completeness of eigenstates 




FIG. 3: Amplitude and density correlators ci(£o;£i) (points) 
and Cff (£o;£i) (error bars) are compared to the semiclassical 
prediction Eq. (|13f> (solid line). An error bar shows mean 
value and standard deviation of those data from Fig.|2i which 
are on the line x\ — xo + po but away from the exceptional 
curves given by Eq. (J2UJ) . 



Finally, this overlap is approximated semiclassically with 
the help of Eqs. ®, Q)- For example we find at t = 1 

M£ )| 2 = 2 3/2 (9 + [2V"(x +p ) ~ 1] 2 )" 1/2 . (13) 

Note that ci(£o;£i) = |&i(£o)| 2 depends only on x\ = 
xo + po which explains the stripes in Fig. 03. In Fig. |3| 
we check the accuracy of Eq. (|T3]) by comparing it to the 
exact value for N = 1024 (points vs solid line). 

We continue with a semiclassical theory for the den- 
sity correlator Eq. ((10)) which allows to understand the 
main features observed in Figs. [2]3,d. For this purpose 
we rephrase Eq. (|12|) as 



U'Wo) 



: a t \a t ) +r t \pt) 



(14) 



where \p t ) is a normalized state which is localized near 
£t but orthogonal to the coherent state at this position, 
(pt\ot.t) = 0- In other words, \ctt) and \pt) can be con- 
sidered as part of some orthonormal basis spanning the 
Hilbert space. In the spirit of random-matrix theory 
one can assume that the coefficients of an eigenstate 
\n) in this basis, and in particular a ty7l = (n\a t ) and 

Pt.n — 

(n\p t ) are uncorrelated, 



P(pt ) = P(pt <n ) P(a t>n ) ■ 



(15) 



Deferring a discussion of the validity of Eq. ((15)) we first 
point out its implications. We multiply the identity 



Ho 



KnlU^ao)? 



• 2Re a* t r t a* t n p t , n + \r t \ 2 \pt,-, 



(16) 

|2 



with H t ^ n = |a^ n | 2 and average over n. Then the second 
term on the r.h.s. vanishes since, according to Eq. ((15)) . 
the phases from a\ n and p t , n are uncorrelated. Also the 
we find £) n e +1£nt/n h^ (£o)hn(0 = (a^|C/*|a^ ). As the third term factorizes and gives (|o^ n | 2 ) n |n| 2 (|p^ n | 2 ) n = 
semiclassical limit is approached, U^a^) localizes at f t (1 - \a t \ 2 )/N 2 after using the normalization implied by 
and the overlap with |c^) vanishes for any other given Eq. JTIJ). Substitution of these results into Eq. jnj yields 
point (Fig. 2a). Therefore all relevant correlations are after straightforward calculation 
given by c t (f ;6) = M£o)| 2 with 



(12) 



C H fa',Zt) = \a t \ 2 V(NI t -l)/(Nlo-l) (17) 
~ M 2 , (18) 
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where X t = J2 n ^tn 18 the inverse participation number 
in Husimi representation. The approximation Zq ~ T t 
leading to from ([T7)) to ([T%|) is justified since deviations 
from the constant X = 2/N expected within RMT are 
mainly due to the influence of periodic orbits |3( which 
is approximately equal if two points £o and £ t connected 
by a short trajectory. 

Eq. (|T%|) suggests that also C# (£o ; £ 1 ) * s given by 
Eq. (p"3j) . This is confirmed by the data in Figs.[2|i and|2l 
(error bars vs solid line) if the conspicuous sharp curves 
of very large correlation (Fig. 01) are ignored. To un- 
derstand the origin of such exceptional points where the 
above theory fails we return now to the crucial step in 
the derivation, Eq. (|15j) . A RMT assumption like this is 
justified only after all important non-generic effects have 
been accounted for. These include in particular symme- 
tries and semiclassical contributions from short trajecto- 
ries. Hence, as \a t ) and \p t ) are both concentrated near 
£t, we expect that Eq. (|T5|) breaks down whenever this 
point is on a symmetry line or on a periodic orbit. For 
example, if £o = £t we have obviously C#(£o;6) = 1 
but | at | 2 <C 1 if the orbit is sufficiently unstable. In our 
model the only symmetry is time reversal: the unitary 
transformation U = D^UD with D = exp(— iV/2h) leads 
to a symmetric propagator U = U T . As a consequence 
the eigenfunctions of U are real in position representa- 
tion and their Husimi density is symmetric with respect 
to p = 0. After back transformation to the original basis 
we infer from this symmetry that the Husimi density at 
£ is correlated with the density at 

i=(x,-p-V'(x)) (19) 

which is confirmed in Fig. 03. For the torus map Eq. (0) 
the point £ t is invariant under this transformation if [2p t + 
V'(x t )] mod 1 = 0. Indeed we obtain from this condition 
with t = 1 the correct functional equation for the lines of 
excess correlation in Fig. 01, 

[2p - V'(z +Po)]modl = 0. (20) 

Note that short periodic orbits of period 1 and 2 must 
be invariant under time reversal and are located on these 
lines. In other numerical data for maps without time- 
reversal invariance we found excess correlations only at 
isolated points corresponding to the short orbits. A quan- 
titative and more detailed analysis of excess correlations 
is still to be done. We expect that this will be feasible by 
applying Gaussian wavepacket dynamics again, now in 
order to correlate the states \a t ) and \p t ) semiclassically. 

We end this paper with a remark on the parameter 
| at | 2 which essentially determines the Husimi correla- 
tions. Fig. 03 shows that it can be very large even in a 
completely chaotic system. In fact |ai| 2 (Eq. (|13|)) is not 
directly related to the local stability eigenvalues of the 



classical map. As explained after Eq. @, this changes 
only when the width of the Gaussian is much larger than 
that of a coherent state, i.e., when \a t \ 2 is already very 
small. 
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